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Abstract 

We analyze the refocusing properties of time reversed waves that propagate in two different 
media during the forward and backward stages of a time-reversal experiment. We consider two 
regimes of wave propagation modeled by the paraxial wave equation with a smooth random re- 
fraction coefficient and the Ito-Schrodingcr equation, respectively. In both regimes, we rigorously 
characterize the refocused signal in the high frequency limit and show that it is statistically 
stable, that is, independent of the realizations of the two media. The analysis is based on a 
characterization of the high frequency limit of the Wigner transform of two fields propagating in 
different media. 

The refocusing quality of the back-propagated signal is determined by the cross correlation of 
the two media. When the two media decorrelate, two distinct de-focusing effects are observed. 
The first one is a purely absorbing effect due to the loss of coherence at a fixed frequency. The 
second one is a phase modulation effect of the refocused signal at each frequency. This causes 
de-focusing of the back-propagated signal in the time domain. 

1 Introduction 

The refocusing of back-propagated pulses in time-reversal experiments has attracted a lot of attention 
recently both in the physical and mathematical literature; see [UEl El E3 EEU ECU Ell E3 and references 
therein. A time reversal experiment consists of two stages. In the first stage, a signal is sent from 
a localized source to an array of receiver-transducers that record the signal in time. In the second 
stage, the signal is time reversed and re-emitted into the medium, that is, the part that is recorded 
first is sent back last and vice versa. It has been observed experimentally and justified theoretically 
that the back-propagated signal refocuses much more tightly at the location of the original source 
when propagation occurs in a highly heterogeneous medium rather than in a homogeneous medium. 
Moreover, the shape of the back-propagated signal does not depend, under appropriate assumptions, 
on the realization of the underlying medium if it is modeled as a random medium. 

In order to obtain a tight refocusing, it is important that the underlying media do not change 
during the two stages of the time reversal experiment. Several experimental studies have demon- 
strated that the refocusing of time reversed waves degrades as the back-propagating medium is 
modified |19| I25j. The modifications in the refocusing properties have been analyzed in (Hj in the 
weak coupling regime based on the formal theories of radiative transfer and diffusion equations for 
time reversed waves propagating in random media They have also been rigorously analyzed in 
the one-dimensional setting in the regime of strong fluctuations and wave localization. It has been 

'Department of Applied Physics and Applied Mathematics, Columbia University, New York NY, 10027; 
gb2030@columbia.edu 

^Department of Mathematics, University of Chicago, Chicago IL, 60637; ryzhik@math.uchicago.edu 



shown in £Q that the re-propagated signal is both not as tightly focused and is no longer statistically 
stable when the two media are different in the one-dimensional case. 

Here we consider time reversal in changing media for two models of wave propagation: the 
paraxial regime and its white noise limit. These regimes model multi-dimensional propagation of 
wave pulses with beam-like structure so that backscattering in the main direction of propagation of 
the beam can be neglected. Time reversed waves in these regimes have been analyzed in pfl IB1 12*2*]. 
We characterize the modifications incurred in the radiative transfer equations modeling time reversal 
as the medium of back-propagation changes. They are described in terms of the cross-correlation of 
the two media of propagation and are similar to those derived formally in We also show that 
the back-propagated signal is still statistically stable, that is, independent of the realizations of the 
random media provided that the correlation functions remain the same. This is similar to what 
was obtained in I22| in the case when the two media are identical and is consistent with the 
numerical simulations in [HJ. This contrasts, however, with the results obtained in the localization 
regime in pQ, where statistical instability has been demonstrated in one dimension. 

As in the pioneering paper on the multi-dimensional time reversal |12j . the characterization of 
the back-propagated signal in the high frequency limit is carried out by analyzing the correlation 
function and the Wigner transform of two wave fields. The main novelty is that we now consider the 
Wigner transform of two fields propagating in two different media |181 123] . Time reversal is the first 
application where such correlations seem to be of a practical interest. Our theoretical analysis is very 
similar to that in [§] and is based on the construction of approximate martingales and perturbed 
test functions. 

The rest of the paper is organized as follows. Section [2] presents the equations modeling time 
reversal in changing media in the paraxial regime. The main results on the characterization of the 
time reversed signal in the paraxial regime are given in Section |31 The theory in the Ito-Schrodinger 
regime is carried out in Section^ In both cases, we observe that the focusing of the back-propagated 
signal at the original source location deteriorates as the cross-correlation of the two media decreases. 
This de-correlation is analyzed in detail in Section Section El offers some concluding remarks. 

Acknowledgment. This work was supported by ONR grant N00014-02- 1-0089, DARPA-ONR 
grant N00014-04-1-0224, NSF Grants DMS-0239097 (GB) and DMS-0203537 (LR), and two Alfred 
P. Sloan Fellowships. 

2 Two-media Time reversal in the paraxial regime 

In this section we generalize the time reversal setting presented in [@] to the situation where the 
media differ during the forward and backward propagation stages. 

2.1 Paraxial wave equation and scaling 

Propagation of acoustic waves is described by the scalar wave equation for the pressure field p(z, x, t) 



Here c(z, x) is the local wave speed, which we model as a random process, and the Laplacian A is 
both in the direction of propagation z and the transverse variable x G M. d . The physical dimension 
is d = 2 although out theory applies to any d > 1. The wave speed c(z,x) is different during the 
forward and backward propagation stages of time reversal. 

The paraxial (or parabolic) approximation of wave propagation consists of assuming that the 
wave field has a "beam-like" structure in the z direction and that back scattering in the z direction 



_L (Pp 

c 2 (z,x) dt 2 



Ap = 0. 
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can be neglected |24l . This implies the approximation 



p(z,x,t)» / e ik( - z - cot ^(z,x,k)codk, (2) 



where the function tp satisfies the Schrodinger equation 

oib 

2ik — (z,x,k) + A x ip(z,x, k) + k 2 (n 2 (z, x) - l)V>(z,x, £;) = 0, ^ 
^(z = 0,x, k) = -0o(x, k) 

and A x is the Laplacian in the variable x. We have defined the refraction index as n(z, x) = cq/c(z, x) 
where cq is a reference speed. Note that (j3J) is an initial value problem in the z-variable. Theoretical 
justifications of the passage from the wave equation to the parabolic approximation can be found in 
□ EH- 

We analyze the high frequency regime, where waves undergo multiple interactions with the 
inhomogeneous medium and wave propagation may be described by macroscopic equations in ap- 
propriate limits. To quantify these limits, we introduce some scaling parameters. Let L x and L z be 
the overall propagation distances. We re-scale x and z as L x x and L z z with the new x and z being 
non-dimensional O(l) quantities. In order for the paraxial approximation Q to be valid one has to 
assume that L x <C L z . 

Let l x and l z , be the transversal and longitudinal correlation lengths of the heterogeneous 
medium. Upon recasting the refraction index as 

n 2 (z,x)-l = -2aV(f,f), (4) 

'•Z ''X 

the above equation © becomes in the re-scaled variables 

2ikdip 1 2 ^ ,L z z L x x 

77 aJ + If ^ ~ 2k aV{ —> 17^ = a (5) 

Let us now assume that the medium and the typical wavelength of the propagating waves satisfy 
the following scaling assumptions: 

£ = f = f«l, kL z = K -(±-\\ a = ^. (6) 

These constraints imply that we are in the high frequency regime when the non-dimensional wave 
number k is of order 0(1). Note that there is one free parameter left in the above relations, namely 

I r = e\ n > o, (?) 

J^z 

where rj > is necessary to be compatible with the paraxial approximation and to ensure that 
L x -C L z . The relations © quantify how the correlation length and the strength of the fluctuations 
are related so that the parabolic wave equation Q in the radiative transfer scaling is given by 

iKe ^L + £: Ax v - k 2 V^v( -, -V = o. (8) 



dz 2 \e e 

The above equation is our model for wave propagation in this section. We will see a different 
scaling in Section |1J This equation is a Schrodinger equation with "time" -dependent potential, as 
the potential depends here also on the variable z. 
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The above choice of scaling implies that 

T Z =T Z =£V<<1 > (9) 

so that the medium is physically anisotropic: fluctuations in the longitudinal and transversal direc- 
tions are not defined at the same scale. Only in the limit L x /L z — > 1, i.e., r/ — > do we recover 
a statistically isotropic medium. This limit, which is more relevant in many practical problems, is 
much more difficult to handle mathematically EH] . The paraxial approximation in the radiative 
transfer regime presented in this section shares most of the physical aspects of the isotropic model 
and is much more amenable to a rigorous mathematical treatment. 



2.2 Time reversal modeling 

We are interested in the refocusing properties of tightly localized pulses. We assume that the center 
of our pulse is a point xo and that its spatial width is e, so that the typical wavelength in the system 
is e. We thus scale our initial condition for the Schrddinger equation as 

V>(s = 0,x,«) = z/> (^^,A (10) 



During the forward propagation phase, we assume that the medium is described by fluctuations 
Vi(z,x). The Green function associated to (jBJ) is then the unique solution to 

. <9G/-(z,x, K;y) e 2 9 (z x\ 

ike — h yA x G/(z,x, K\y) - K^eVil^-, -JGf{z,x,K;y) = q-q 

G/(0,x, rc;y) = 6(x - y). 

Let us assume that waves propagate for a distance z = L = cqT along the z axis, or equivalently for 
a time T. The solution at z = L is given by 

^_(L,x,k) = / Gf(L, x, k; y)ipo ( — — —,Kjdy. (12) 

The signal is then recorded on a domain of small (but of order O(l)) aperture - this is modeled 
by multiplication of the signal by a compactly supported function x( x )- We also allow for some 
blurring at the detectors so that the re-emitted signal after time reversal is given by 

V> + (L,x,«)=x(x) f £ - d f(^^) X (yW-(L,y,K)dy. (13) 



Here * denotes complex conjugation and corresponds to time reversal. Indeed, the time reversal 
t — > — t in the time domain amounts to complex conjugation e lu)t — > e~ lujt in the frequency domain. 
The blurring must be controlled at the scale of the wavelength e for otherwise all the coherent signal 
would be irretrievably lost. 

It now remains to model back-propagation to the hyperplane z = 0, that is, again during a 
time T. The back-propagation takes place in a different medium described by the random potential 
Viiz-, x) whose Green's function satisfies 

dGb(z, x, k; y) e 2 , . 9 /z x\ . 

we — h — A x G 6 (z, x, k; y) - K A ^eV 2 y-, -J G b {z, x, k; y) = ^ 

G fe (0,x,«;y) = 8(x - y). 
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After back-propagation for a distance L along the z-axis and a second time reversion (complex 
conjugation in the frequency domain) we obtain that the re-propagated signal takes the form 

$ b (x,k)= [ G* b (L,x, K ;y)iP*+(L,y,K)dy. (15) 

JR d 

The second conjugation is performed so that when full measurements are available, that is, % = 1, 
and the detectors are perfect, so that /(x) = <5(x), we recover the original signal exactly: V ;B ( X ) K ) = 
ip(z = 0, x, k). 

We are interested in the back-propagated signal in the vicinity of xo and define 

^f(€,K;xo) =$ B (xo + e£,K). (16) 

Summarizing the successive steps described above, we can relate the back-propagated signal to the 
initial signal as 

^f(£,«;xo)= / Gl(L,y^ + e£,K]r])Gf(L,yti + eC,K,y)x(v,y)^o(G,K)dCdridy, (17) 
Jm? d 

where we have used that G(L, x, k; y) = G(L, y, k; x) as can be seen from the equation satisfied by 
the Green function and where we have defined 

X(V,Y) = X(v)x(y)f(—) = xivMOj^ [ f^e^e^^dq. (18) 
\ e J (2ir) a J R d 

The above notation implicitly defines our convention for the Fourier transform /(q) of /(x). We 
observe that the back-propagated signal in 1)17(1 involves the product of two Green's functions at 
nearby points. The Wigner transform is thus a very natural tool to understand the statistical 
properties of this two point correlation Following [51 [7] we introduce the functions <5/,& as 

Q /j6 (L,x,K;q) = / G /j6 (L,x, K;y) X (y) e -^/ £ dy, (19) 

which solve the initial value problems 

f b ^ X\ 

i£K £ (z,x,K;q) + yA x Q /jfe (z,x, K;q) - K 2 ^eV Sf> -j Qf jb (z,x, k; q) = 0, ^ 
Q /ib (/=0,x, K ;q) = X (x)e-^ x / £ . 

We then define the Wigner measure W e as 

W E (z,x,k,K) = / /(q)[/ e (z,x,k,K;q)dq, (21) 

where £/ e is the Wigner transform of the auxiliary functions Qf t b defined by 

U e (z,x,k,K;q) = [ e 4ky Q/(z,x - -jf, K;q L )Qt(z,x+ ^,K;q)-^. (22) 
J M d 2 2 (2ir) a 

The main reason for introducing the above notation is that the back-propagated signal can be recast 
in terms of the Wigner measure as 

^f(€,«;xo)= / e^-^W £ (L, Xo + el±A,k, K )My^)^- (23) 

J R 2d 2 (27T) a 

Thus in order to understand the macroscopic properties of the time reversed signal xftf in the high 
frequency limit, i.e., as e — > 0, it suffices to analyze the Wigner measure W e in the same limit. This 
task is taken up in the following section. 
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3 Stability of waves in changing environment 



3.1 The main result 

We consider in this section the general problem of the correlation of solutions of the linear paraxial 
Schrodinger equations in two different albeit correlated random media. We let tp £ (z,x) and </> e (z,x) 
be the solutions of the family of Cauchy problems 



l£K- 



dips 
dz 



+ -Ail> e -K i V£Vi[-,-)i/>. 



Z X 



£ £ 







(24) 



^(0,x)=^°(x;C) 



and 



z x 



ieK-p- + — A^» e - K 2 y/eV 2 -, - 

OZ 2 \£ £ 

£ (O,x) = 0°(x;C) 



0. 



(25) 



with two different random potentials V\ and V 2 - The initial data depend on an additional random 
variable £ defined over a state space S with a probability measure cfyz(£). It accounts for the 
consideration of a mixture of states rather than the single solution of the Schrodinger equation. The 
mixture of states arises naturally in the time-reversal set-up, because of the integration over the 
wave vector q in 1)21(1 . This introduces additional regularity into the problem, which is crucial to 
obtain statistical stability. 

The cross Wigner transform is defined by 



W e (*,x,k) 



e ik ' y ^ e ( z, x 



i z, x + — ; ( 



(2tt)< 



MO- 



The evolution equation for the Wigner transform is 



OZ K 



3 ip-x/e 



V 1 (- >P )W £ (k 



P )W e k + 



dp 



■ (26) 



Here V(z,p) is the partial Fourier transform of V(z,x.) in x only. We will assume that the initial 



data We(0, x, k) converges strongly in 



to a limit Wq(x, k). This is possible due to 



the introduction of the mixture of states - the integration against the measure [i{d^)- although the 
Wigner transform of a pure state is not uniformly bounded in L 2 (M. d x W 1 ) [20j. The evolution 
equation (|26|) preserves the L 2 -norm so that in order to identify the limit of W e as e — ► 0, it suffices 
to consider initial data 

W £ (0,x,k) = W o (x,k) (27) 

that are independent of the parameter e. 

We assume that the random processes V\p{z) are statistically homogeneous in space x and "time" 
z, have mean zero and rapidly decaying correlation functions Rij(s,y): 



E{V t (z,*)} = 0, E{V i (z + s,-K + y)V j (z,x)} = R ij (s,y), 
We denote by Rij(u,p) the corresponding power spectra: 



i,j = 1,2. 



E 



{^(u;,p)^V,q)} = (27r) d+1 ^( W ,p)5( W + W ')5(p + q), ^>,p)= / e^^R^t, x)dtdx. 
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We will also assume that the partial Fourier transforms Vj(z, p) in x only are almost surely supported 
in a deterministic compact set { ||p|| < C} and the total mass is also almost surely uniformly bounded: 



J \dVj(z,p)\ < C, 



with a deterministic constant C. We denote the state space of such spectral measure by V. 

We further assume that the joint random process V(z) = (Vi(z), ^(z)) is Markovian in the 
variable z with a generator Q (written in the Fourier domain) that is bounded on L°°(V), has a 
unique invariant measure tt(V) and a spectral gap a > 0. This means that 

Q*tt = 0, 

and if (g, tt) = 0, then 

\\e rQ g\\ L ™ < C\\g\\ L ™e- ar . (28) 
Given (|28[) . the Fredholm alternative holds for the Poisson equation 

Qf = 9, 

provided that g satisfies (ir,g) =0. It has a unique solution / with (tt, f) = and < C||fl'||i^ > - 

The solution / is given explicitly by 

POD 

f(V) = - / dre^g(V), 
Jo 

and the integral converges absolutely because of ((2*5|) . 

The main result of this section is that under the above assumptions, the following theorem holds. 
Let us define the operator 



£/(x,k) 



a , P 2 -k 2 , w ( v ^ii(^, P -k)+i? 22 (^,p-k) w 
^i2( 2 ' p ~ k )^o(p) 2 Wb(k) 



dp 



with 



(2ir) d 

-OT(k)Wo(k) (29) 

mn 1 [ j [ d P ^22(r-,p) - #ii(r,p) r . . . . , . . 
II(k) = - dr ——T exp{ir(k - p/2) • p)sgn(r) 

1 JR JRd ( 27r ) 2 

f f #22(w,k- p) - #n(w,k - p) cL;dp 

= / p - v - / | D |2_i k |2 ^Tw+i • ( 30 ) 

Jr* Jr co- |p| 2 |k| (2vr) d+i 

Here, i?(r, p) is the partial Fourier transform of R in x only. We denote the standard inner product 
on L 2 (M? d ) by (/, g) = f R2 d /(x, k)p(x, k)dxcfk. Then we have the following result. 

Theorem 3.1 Under the above assumptions, the Wigner distribution W e converges in probability 
and weakly in L 2 (M. 2d ) to the solution W of the transport equation 

dW n 

k— + k • V^W = n 2 CW. (31) 

oz 

More precisely, for any test function A £ L 2 (M. 2d ) the process (W s (z),X) converges to (W(z),X) in 
probability as e — > 0, uniformly on finite intervals < z < Z. 



7 



3.2 Proof of Theorem EO 

The strategy of the proof is very similar to that in [§]. Observe first that since the Wigner equation 
preserves the L 2 -norm, the joint process (W £ (z),V(z)) is a Markov process on X x V, where X = 
{||W||2 — C} is an appropriate ball in L 2 (R d x M ). The corresponding family of measures P £ on 
the right-continuous paths on X is tight, as can be shown in a way identical to |5j and |S] (see also 
jllj for a detailed calculation in a similar setting). 

Given a test function A(z,x, k) we will show that the functional 

G\(z) = (W, A) - jT /w, f ^ + • V x + A^ (32) 

is an approximate P e -martingale. More precisely, we show that 

|E ft {G x [W](z)\F s } - G x [W](s)\ < C x ,zV^ (33) 

uniformly for all W G C([0,Z];X) and < s < z < Z, with a deterministic constant C,\z- The 
weak convergence of the probability measures P e together with l)33|) imply that E{VF e } converges to 
W. In order to establish 1)33(1 we will construct another functional G\ that is an exact martingale 
and that is uniformly close to G\. This is done by the perturbed test function method. A similar 
argument applied to (W, A) 2 implies that E{VF £ <8> W £ } converges weakly to W <8) W. This implies 
convergence in probability. In order to simplify the notation we set k = 1 throughout the proof. 

Step 1. Convergence of the expectation. Given a function F(W,V) let us define the 
conditional expectation 



E Pe 



W,V ,z 



{f(W, V)} (r) = E Pe |f(W(t), V"(t))| W{z) = W, V{z) = , r > z, 



where P e is the joint probability measure of V and W e . The weak form of the infinitesimal generator 
of the Markov process generated by V\ 2 and W £ is given by 



d 

hence 



KyA {WA{V)) } {Z + k) 



h=0 



±{W, QX) + (V, + k • V x - \k.[V, *] ) 



G e A = (W, \{V)){z) - J* (w, Q-Q + A + k • V x - i=/C[y, *]J A^ (s)ds (35) 
is a martingale. The skew-symmetric operator K, is defined by 

!C[V, m ^KV) = - [ §Me^(x,^,k-f)-i/ ^P) e ^( x ^,k + f). (36) 

The generator (|34|) results from the Wigner equation written in the form 

8W 1 ~ z x 

+ k • V X W £ = -=K\y(-), -}W £ . (37) 
oz y/e e e 

The following lemma is the key element to show that E{W e } — > W, solution of (|31|). 

Lemma 3.2 Let X(z, x, k) G C 1 ([0,Z];5) 6e a deterministic test function, and let the junctionals 
G\ and G\ be defined by \3ty) and \35)) . respectively. There exists a deterministic constant C\ > 
and a family of perturbed random test functions \ £ so that ||A £ — A 1 1 2 < C\^fe almost surely and 

\\G\ E {z)-G x {z)\\ Lao{v) <C x V~e (38) 
uniformly for all distances z G [0, Z] . 
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The proof of this lemma is presented in Appendix The weak convergence of the probability 
measures P e and Lemma 13.21 imply that E{W e } — > W, weak solution of 



{W(z),X(z))-{W ,X(0)) 



ds( W 



0_ 

ds 



+ k- V x + £* J A) (a) = 



(39) 



which is nothing but the weak form of (|31j>. 

Step 2. Convergence in probability. We now look at the second moment E{(W e ,A) 2 } 
and show that it converges to (W, A) 2 . This implies convergence in probability. The calculation 
is similar to that for E{(W £ ,A)} and is based on constructing an approximate martingale for the 
functional (W ® W, [/,), where /i(2,xi,ki,X2,k2) is a test function, and W ® W(z, xi,ki,X2,k2) = 
M / (z,xi,ki)W / (z,X2,k2). As before we consider functionals of W and V of the form F(W,V) = 
(W ® W,[i(V)), where /i. is a given function. The infinitesimal generator acts on such functions as 



dh w,v,z 



(W®W,n{V))} (z + h) 



h=0 



-(W<g> W,Q\) + (W®W,H £ 2 fi), 



where 



with 



and 



j=l v 



V, 



x' 



/'• 



(40) 



(41) 



* jR d ( 27r ) 2 



1 /■ dV r 2 (p) ) 



! < ^V(ki + ?, k 2 



Therefore the functional 



^(P) c i(p-^) urk| k , _ £ 



i jRd (2vr) d 2 



2,s 



(iy® w,n<y)){z) 

W®W, (-Q + — + ki • V X1 + k 2 • V X2 - -^(Ki[V, — 
\e oz yJE e " 

is a P e martingale. The following lemma is proved in Appendix iBl 



(42) 



x 2 



Lemma 3.3 Let fi(z, x 1; k l5 x 2 , k 2 ) be a deterministic test function and let the functional G]f be 
defined by \4^D - Then there exists a deterministic constant C > so that 



\G 



2,s 



G 2 /\ < CV~e 



with 



G 2 f = (W®W,n)(z) 



w ® w > iL + ki ■ Vxi + k2 ■ Vxa + 



/x ) (s)ds 



(43) 



(44) 



and mi/t a deterministic operator £ 2j£ swc/i that ||£ 2 e — £* ® £*||x2^^ 2 — * as e ^ 0. 

Lemma f3 . 3l implies immediately that for any test function [i we have E {(W E ® W e ,fj)} — > {W<S>W, jj). 
If we take /i = A ® A we get E{(W £) A)} -> (TF, A) 2 and hence (W e , A) -> (W, A) in probability. This 
finishes the proof of Theorem 13.11 
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4 The Ito-Schrodinger regime 



We consider in this section the regime where the ratio l z /L z of the correlation length l z of the 
fluctuations in the z direction to the propagation distance L z is the smallest parameter in the 
system. 

4.1 Ito-Schrodinger equation 

Let us recall the Schrodinger equation © 

2ik dib 1 . 9 ,L z z L r x x , 

The scaling assumptions © are now replaced by 

e = — < 1, — = e^ a , a > 0, kL z = -[ — J , o- = e 2 _. (46) 

The constraint a > indeed implies that l z /L z is smaller than any other dimensionless term in the 
system. With these assumptions, (|45|) may be recast as 

Because the variations in z of the potential are faster than any other quantity in the above equation, 
we can formally replace 

^V(^7) dz ^ inB(dz^), (48) 
its white noise limit, where B(dz,x) is the Wiener measure described by the statistics 

E{£(x, z)B(y, z')} = K(x - y)z A z'. (49) 

Here, E{-} means mathematical expectation with respect to the Wiener measure, K(x) is the corre- 
lation function of the random fluctuations and zAz' = min(z, z'). The paraxial Schrodinger equation 
then becomes the following stochastic equation 

IS X! 

dip(z, x) = — A x ^(z, x.)dz + iKtp(z, x) o B(dz, —). (50) 

Here, the notation o means that the stochastic equation is understood in the Stratonovich sense 
|17l 121) . In the ltd formalism, it becomes the following Ito-Schrodinger equation 

1 fie 2 T w„x\ . , / x , . . , / n , x, 



di/j(z, x) = k 2 K(0) A x ip(z, x)dz + inip(z, x)B(dz, -). (51) 

2 V K J £ 

We do not justify the derivation of (|51|) here. It was shown in [2] that the paraxial approximation 
and the white noise limit can be taken consistently in the one-dimensional case. 

As in the paraxial regime, we still have one parameter left, namely L x /L z , which we choose as 
in lf?jl. We then verify that 

l_x_ _ jx_ Lx_ Az _ £ T}-a /^2) 
l z L x L z lz 

Thus with the choice rj = a, the Ito-Schrodinger equation ((5T|) can be used to model isotropic 
fluctuations. 
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4.2 Time reversed waves in changing media 

The formalism presented in Section 12.21 applies in the white noise limit as well. We can still define 
the functions Q/,&, which now solve 



dQ/ )6 (z,x,K;q) = -y— - K 2 K li2 (0)jA x Qf, b (z,x,K;q)dz + inQf^z, x, re; q)B l)2 (dz, -), ^ 
Q/, 6 (0,x,re;q) =x(x)e- iX ' ^ £ , 

where the Wiener measures B\ )2 are described by different statistics K\^ 2 for the forward propagation 
(index 1) and the backward propagation (index 2). The cross- correlation of the two media, is defined 
by 

E{B m (x, z)B n (y, z')} = K mn (x - y)z A z', 1 < m, 7i < 2. (54) 

We will see in what follows that the relative strength of the cross-correlation K\ 2 compared to the 
auto-correlation functions K mm determines the quality of time-reversal. 
Upon defining 

U e (z,x,k, re;q) = / e' ky Q / (z, x - k; q)Q^(z, x + k; q)-^, (55) 
J Rd J 2 2 (2vr) a 

as in (1221) and 



W e (z,x,k,re) = / /(q)[/ £ (z,x,k, re;q)dq, (56) 
as in (|21|) . we obtain that the back-propagated signal is given as in Q23JI by 

( (£,re;x ) = / e *«-y) W£ ( L)Xo + e y+l )k) ^ ( y)K) ^. (57) 

J R 2d 2 (27T) a 



V'. 



The high frequency limit of the time reversed signal is thus again modeled by the limit e — > in the 
above equation. 

4.3 High frequency limit of time reversed waves in changing media 

In the high frequency limit, we have the following result 

Theorem 4.1 Let k G K fixed. Let us assume that the initial condition ipo(y, K ) £ -L 2 (M' i ), i/ie /iZier 
/(x) G L 1 (]R ci ) n L 00 ^^), and i/ie recorder function x( x ) «s sufficiently smooth. Then ipf(£, k; xo) 
converges weakly and in probability to the deterministic signal 

V> B (£,re;x ) = / e ik ^(L,x ,k,re)^o(k,re)dk, (58) 
where W(L,xq,\s., re) solves the following radiative transfer equation 



2 JRd 



dW 1 

— + -k- X7 x W + k 

oz re 2 j R d (59) 

TU(0,x,k, K ) = /(k)x 2 (x). 

Moreover for a smooth test function of the form A(£,xo) = A(xo)//(£), we /lave an error estimate of 
the form 

E{(Vf - E{^f}) 2 } < C £ d ||A||| 2(Rd) ||^ |l!2 (Rd) , (60) 
uniformly in L on compact intervals. 
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The main steps of the proof of the theorem are very similar to that in the paraxial regime. However 
the mathematical analysis is substantially simplified by the fact that statistical moments of the field 
ipf and the associated Wigner transform W e satisfy closed-form equations. We refer the reader to 
|14l 1171 I26| for basic results about the stochastic partial differential equation (|57|1 . The proof of the 
above theorem can be carried out as in Pj. We highlight the differences that appear because of the 
change of media during the forward and backward propagation. 
Let ipi and ip 2 satisfy 

dip m (z,x) = -( K 2 ET m (0))A x ^ m (z,x)(i2: + iKip m (z,x)B m (dz, -), to =1,2. (61) 

2 \ k / e 

We define the second moment TO2 (x, y ) as 

m 2 (z,x,y,«) = E{i[i 1 (z,x + ^-,k)iJ)Z(z,x - ^-,«)}- (62) 
By an application of the ltd calculus [U we obtain that 

d(ipi 0, x)^2 {z, y)) = ^i (z, x)dip2 {z, y) + dipi (z, x)ip2 {z, y) + dipi {z, x)dip* 2 (z, y). 

We insert Q61|) into the above formula and taking mathematical expectation, obtain after some 
algebra [3] an equation for m-i'- 



dm 2 



-V x -V y m 2 (z) -k 2 ^ ^— ^ - K 12 (y)jm 2 (z). (63) 



dz 

Now, defining the Wigner transform of the two fields as 



W 12 (z,x,k,K) = -±- s f e ^i(z,x-^, K )^(z,x + ^, K )^, (64) 
(2vr) a J R d 2 2 



we find that 



m 2 (z,x,y,K)= e lky E{W 12 }(z,x,k,K)dk. (65) 
Therefore, E{Wi2} solves the following equation 

^ + ik • V*W + K 2 Kn(0) + K 22 (0) w = k2 I £ i2(p _ k)w{p)d ^ (g6) 

OZ K 2 J K d 

This is the integro-differential equation in Q59|) . By construction, E{C/ e } defined in (|55|). whence 
E{W e } defined in (|56j). satisfy the same equation. 
Let us now consider the fourth-order moment 

TO 4 (z,X,y,Z,t,K) = E{lpi(z,X + ^,«)^(Z,X- ^ptt)'</'l(z,Z + y,K)^2(^,Z ~ ( 67 ) 

We deduce from the application of ltd calculus to four arbitrary functions 

d(ifai/>$ifail>t) = ^ 2 tp3iptdtpi + - ■ ■+ ipitp 2 ip3dipt + ipiip* 2 dip z dipl + - ■ ■ + ip^ldipidip* 2 , 
that rri4 solves the following equation 

—q^- = -(V x - Vy + V$ • Vt)m^{z) -JCm^z), 
£(x,y,£,t) = K n (0) + K 22 (0) - K 12 (y) - K 12 (t) 

11 s/ 2 ' e , 2 y 
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Let us now introduce the second moment of W12' 

W(z,x,p,£,q,re) = W 12 (z,x,p, k)W 12 (z,£,cl, k). (69) 

We verify that 

m A (z, x, y, z, t, k) = / e'^^EjWKz, x, p, y, t, K )dpdq, 

so that E{W} solves the following equation 
dW 1 



(70) 



+ -(p • V x + q • V € )W + k 2 (^h(0) + K 22 (0))W = k z £ 2 W + as=*£i2W, (71) 



where 



£ 2 W = / (^i 2 (p-p , )5(q-q)+^12(p-p)5(q-q , ))W(p / ,q / )dpW 
C 12 W = £e^ u (K 12 (u)(W( P -^q-|)+W(p + |,q+^)) (72) 
-#n(u)W(p - |, q + I ) - K 22 (u)W(p + |, q - |)) du. 

We thus obtain that both 

W £ (>,x,p,£,q, «;k) = E{J7 £ (z,x,p, «;k)J7 e (z,€, q, re;k)} (73) 

where C/ e is defined in (|55|). and 

W 6 (z,x,p,£,q,«) =E{W e (*,x,p,K)W e (*,£,q,K)} (74) 

where W e is defined in (|56|) . satisfy the same radiative transfer equation (|71|) . There is however a 
fundamental difference between the two latter terms, namely that W £ is bounded in L 2 (R 4rf ) at fixed 
n, whereas U £ is not bounded in the same norm at k and k fixed. Indeed, W e (z = 0) is bounded in 
L 2 (R 2d ), which is not the case for U e (z = 0). The results in [31 section 3] show that E{W e }(^) and 
W £ (z), are then bounded in L 2 (M. 2d ) and L 2 (M 4d ) respectively, uniformly in z > 0. More precisely, 
we have 

W £ (0, x, k) = / e- k y/( y ) X (x + ^) X (x - ^)dy. (75) 
JK d ^ z, 

For /(x) and x( x ) sufficiently smooth, Theorem 4.1 of [3] allows us to conclude that 

\\W £ - E{W E (z,x,p,K)}E{W £ (z,t<l,K)}\\ L2(R4d) < Ce d / 2 , (76) 



uniformly on compact sets in z. This comes merely from the observation that L\ 2 defined in (|72|) 
converges to zero as an operator on L 2 . Moreover, (|75|) implies that W e (^ = 0, x, k, k) converges 
strongly to /(k)x 2 (x) as e — ► since x( x ) is smooth. This implies that E{W e (z, x, k, k)} converges 
strongly in L 2 (M? d ) and uniformly in z and K on compact intervals to W(z, X, k, re) solution to (|59|) 
as e — ► (since the L 2 norm is preserved by (|59jl). 

For a test function A 6 L 2 (M 2rf ), the above convergence implies that 

E{((W £ ,X) - (E{W e }, A)) 2 } < C e d / 2 ||A||| 2(R2d) . (77) 

We deduce that (W e ,A) converges in probability to the deterministic number (W, A) as e — ► 0. We 
have thus obtained the (weak) stability of W £ . Then we can pass to the limit e — ► in (|57[) and 
obtain (|58|). This concludes the proof of Theorem 14. II 
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5 Decoherence in time reversal 



The two preceding sections were concerned with the derivation of the radiative transfer equations 
modeling time reversal when the medium during the backward propagation phase differs from the 
medium during the forward propagation stage. In both regimes we observe that the main quantity 
governing refocusing is the ratio of the cross-correlation terms R\2 and K12 to the auto-correlations 
Rmm and K mm , m = 1,2. When that ratio is large, time reversal refocusing works as if both 
media were the same. When the cross-correlation is small, the coherent effects that produce strong 
refocusing in time reversal are no longer present. 

Let us focus on the two-media effect in the Ito-Schrodinger regime first. We recast ()58|) in the 
Fourier domain and obtain 

^ B (k, k; x ) = W(L, x , k, «)^,(k, k). (78) 

Therefore, the medium acts as a filter between the original signal ^o(k, k) and the refocused signal 
ip B (k, k; xo). The back-propagated signal is all the tighter around xo that the filter is close to a 
constant (in k) non-zero value. Since W satisfies a radiative transfer equation, the regularity of W 
is improved due to the scattering term on the right hand side in (|59[). as is discussed in detail in 
[7j. Indeed, multiple scattering has a regularizing effect. As the change in the propagating media 
increases, the cross correlation K\2 decreases. This weakens the scattering term in (|59)1 . hence 
diminishes the regularizing effect of (|59|) and the re-focusing properties of the time reversed signal. 
Let us assume that K12 is real-valued to simplify the presentation. The weakened refocusing can be 
quantified by recasting the radiative transfer equation (|59|) as 

?K + I k . VxW + K 2 aaW = K 2 f K 12 (p - k)(W(p) - W(k))dp 
dz k J R d (79) 

W(0,x,k,«) =/>) X 2 (x), 
where we have defined the apparent absorption coefficient 

^(p-k M p. ,80, 



As the media decorrelate, the absorption coefficient a a increases up to the value ±(ifn(0) +#22(0)) 
when the two media become completely uncorrelated. The right-hand side in (|79|) then vanishes and 
the back-propagated signal is the poorly refocused signal one would obtain in a homogeneous medium 
with constant wave speed c = cq, albeit with a decreased amplitude by a factor e~ K L ( K ^(°)+ K ^W) . 

Similarly, a signal that is back-propagated in a homogeneous medium would be modeled by 
V2 = 0, which implies that K12 = K22 = 0. So the back-propagated signal would similarly be, up to 
a factor e~ K LKll \°\ the poorly refocused signal one would obtain in a homogeneous medium. Unless 
we have a sufficiently accurate knowledge of the underlying medium, back-propagating a recorded 
signal in a homogeneous medium, for instance on a computer, will not tightly refocus at the original 
location of the source term. 

The situation is somewhat richer in the paraxial regime. The radiative transfer equation takes 
the form 

f)W 1 r ~ n 2 — k 2 

_ + _ k • V X W + K 2 (a a (k) + Hl(k))W = K 2 J d Ri2(^— , P - k)(W(p) - W(k))dp 
W(0,x,k,«) = /(k) X 2 (x), 
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where we have defined 

2 i 2 2 l 2 2 i 2 

^a(k) = y ji(ij n (P-^— ,p-k)+£ 22 (5-^_ ,p- k ))-fi 12 (5-^_ ,p-k)]dp, 

fl22(w,k-p)-flu(w,k-p) ciwdp ( 82 ) 



n(k) = / p.v. / 



IpP-I^I 2 (27r) d+1 



Still assuming that i?i 2 is real- valued, we obtain that <r a (k) is an apparent non- negative absorption 
coefficient and iIT(k) is a purely imaginary modulation term. 

We have seen the role of the absorption a a in the Ito-Schrodinger regime. The role of the new 
modulation term iIT(k) is somewhat different. It also reduces the strength of the right hand side 
in (|81|) but only in the time domain, when we integrate over all frequencies. Let us assume for 
instance that Il(k) is constant. We then verify that W(z) = e tK nz U(z), where U(z) satisfies the 
same equation (|8Tj) with IT replaced by zero. Consequently, II has a tendency to modulate the filter 
W(z) that appears in (|78j) . The modulation is independent of the wave vector k or the position Xo- 
However, it depends on the longitudinal length z and on the reduced wave number k. Therefore, in 
the time dependent time reversal experiments, where the refocused signal p s (0,x, t) is given by (|2*Jl 
with tp replaced by ip B , that is, as an average over reduced wave numbers k (after an appropriate 
re-scaling), the modulation factor II will imply that the back-propagated signal is given by 

p B (0,tt)~ [ e - mco ^ B (0,£,K) CoC k= f e- iKCot e iK2uL Wo(L,x ,k,K)4, (KK)codK, (83) 
it it 

where Wq is the filter obtained when II = 0. Obviously, the magnitude of the above oscillatory 
integral decreases as II increases. The interpretation of the modulation term IT is thus the following. 
Although it does not modify the intensity of the filter W{L, x, k, k) at a fixed frequency, it introduces 
a modulation of order e lK UL that significantly reduces the back-propagated signal recorded in the 
time domain. 

Let us conclude with a remark on the comparison between the radiative transfer equations in 
the paraxial and Ito-Schrodinger regimes. The latter regime should be seen as a limit of the former 
as the oscillations in the z direction become faster and faster. Indeed, the fast oscillations in the 
variable z imply a decorrelation in the term R(x,z), which converges to K(k)5(z). This in turn is 
consistent with R(u>,p) converging to K(p)- It remains to observe that the Hilbert transform (the 
principal value integral in (|82j0 of a constant function vanishes to conclude that IT(k) vanishes in 
the limit of fast oscillations in the z direction. This implies that the oscillatory integral obtained in 
(|8H|) can only be observed in media where the oscillations in the z variable have a sufficiently large 
correlation length. 

All the effects mentioned in this section are in agreement with the radiative transfer and diffusion 
numerical simulations performed in jH] in the so-called weak-coupling regime, which is the limit 
L x w L z of the two regimes considered in this paper and for which no rigorous mathematical 
derivation is available. 



6 Conclusions 

When the medium is fixed during the forward and backward stages of a time reversal experiment, 
the refocusing of the back-propagated pulse is characterized in many high frequency regimes by 
a radiative transfer equation. The solution to the radiative transfer equation acts as a transfer 
function and indicates how the shape of the original source term is modified by the time reversal 
experiment. We have shown in this paper that this picture remains valid when the two media during 



15 



the forward and backward stages differ. We have also described how the constitutive parameters of 
the radiative transfer equation change as the back-propagation medium is modified. Moreover, these 
parameters only depend on the correlation function of the two media. Finally, we have observed 
that the refocused signal was essentially independent of the realization of the random medium. More 
precisely we have shown that the back-propagated signal converges weakly and in probability to a 
deterministic function in the high frequency limit. This results from a similar convergence property 
for the properly regularized Wigner transform of two fields propagating in two different media. 

As the two media are increasingly decorrelated, the refocusing of the back-propagated pulse 
degrades. Two mechanisms are responsible for this degradation. The first mechanism consists of a 
purely absorbing term indicating that wave mixing by scattering is less efficient as the two media 
become less correlated. This effect, though frequency-dependent, can be observed at all frequencies, 
hence also in the time domain. The second mechanism, which is absent in the Ito-Schrodinger 
regime, is a phase modulation phenomenon in the frequency domain. The signal at frequency c^k 
is modified by a phase proportional to k 2 , which has an important cancellation effect in the time 
domain after Fourier transforms are performed. 



A The proof of Lemma 13.21 

Given a test function A(z,x, k) 6 C 1 ([0,Z];5) we define the following approximation 

X E (z, x, k, V) = X{z, x, k) + y/eXl(z, x, k, V ) + eX £ 2 {z, x, k, V) (84) 

with Af 2 C- 2 ') bounded in L°°(V; L 2 (M. 2d )) uniformly in z £ [0, Z\. The functions Af 2 W1 U be chosen 
in such a way that 

\\GUz)-G x {z)\\ LOO{v) <C x V~e (85) 

for all times z £ [0, Z\. Here the functional G £ is defined by (|35|) and the functional G by ()32j) . 

The functions Af and A| are constructed as follows. Let Xi(z, x, k, V) be the mean-zero solution 
of the Poisson equation 

k • V^Ai + QAi = K.X. (86) 

It is given explicitly by 

X 1 (z,x,^,k,V) = - 1 - ^ 'dre* [ ^- e ^^X(z,x,k-^) (87) 

* Jo JR d ( 27r r 2 

Then we let Xf(z, x, k, V") = X±(z, x, x/e, k, V). Furthermore, the second order corrector is given by 
A|(z, x, k, V) = X2(z, x, x/e, k, V) where X2&, x, £, k, V) is the mean-zero solution of 

k- V 4 A 2 + QA 2 = Oi -E{Oi}. (88) 

A mean-zero solution of (|88jl exists according to the Fredholm alternative, as the operator Q has a 
spectral gap. A straightforward calculation presented below shows that 

E{Oi} = -£*A. (89) 

Hence the second corrector is given by 

f 00 r i 

A 2 (z,x,£,k,y) = - / dre rQ £*X(z, x, k) + [Oj(z, x, ^ + rk, k, V) . 

Jo 1 J 
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The above computation and straightforward estimates, as in (Hj, show that 



-A 



where ||Cg H2 < C, with a deterministic constant C > 0. It follows that given by 



GUt) = (W(t),Xe 



ds(\V, (J- s +k-V x + jC*) a\ (*)->/£ £ds(W,g){s) (90) 



is a martingale with respect to the measure -P e defined on D([0, Z]; X x V), the space of right- 
continuous paths with left-side limits. In order to show that (|89|) holds let us compute 

«™-«{-7jC^^«- k -!) + TjL^^«- k + ?)} 

= h + 1 2 + II\ + Ih- 
We compute the four terms above separately: 



I = E 



with 



h = -E\ [ f^e^< r dre* [ ^M e ir((k-p/a).q)+i(€.q) A( k _ £ _ 

y Rd (27t)^ y V W 2 2' 

= ~ So dT I ^i( r 'P) e "' ir((k " P/2)P)A (^ x ' k )(|^J 

^p TT i?n(w,p)A(z,x,k) drexp{ir[w-(k-p/2)-p]}. 



The second term is 



h = E i I ^4e^ 



(27T) 



</ 



/ ^^ e -((k-p/2). q)+ ^. q)A(Z)X)k _P + q ; 
Jo JK d (27r) 2 2 



= ^°°dr I ^ 12 ( r ,p) e —((k- P /2).p) A(Z)X)k _ p) _^_ 

/ (2^+! ^ 12 ^' P) A ( Z ' x ' k ~ P) J q drexp{ir[u- (k-p/2) • p]}. 



The term /I is given by 
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with 



and 



III — E < f ^PljP< r dre rQ I >((k+p/2)-q)+»(£-q)w k , P _ 

= f°dr / i? 21 (r,p)e— « k +P/ 2 )P)A(t,x,k + p)-^ 
J0 JR d ( 27r ) 

(27r) d+1 ^ 21 ^' P) A ( Z ' x ' k + P) y o dresxp{ir[u- (k + p/2) • p]} 
(2 7r )d+i ^ 12 ^' p ) A ( z ' x ' k - P) J dr exp{ir[-w + (k - p/2) • p]} 
(27r) d+1 "^ 12 ^' P ^^' x ' k ~ p ) _/ dr ex P{ ir I w - ( k - P/ 2 ) • P]} 

//2 = -E / / ^M e ^ dre^ / ^M e -((k + p/2)-q) +i (^ q ) A( k + P + H> 
\J Rd (2vr)^ y 7 Rd (2tt)^ V 2 2' 

= - f°dr / i? 22 (r,p)e-« k+ P/ 2 )-P)A(z,x,k)-^ J 

= ~/ (2^)rf+i ^2(^, p)A(z, x, k) J drexp{ir[uj - (k + p/2) -p]}. 



Observe that 



/d d P ®q 

^ 2 —p TT i?i 2 (cu,p)A(2;,x,k - p) y drexp{ir[w - (k - p/2) • p]} 

/dp /" - k 2 — p 2 dp 
Ru({k - p/2) • p, p)X(z, x, k - P)^p = y #12( ^ , k - p)A(z, x, p N 



(27T) 



Furthermore, we also have 



-[/ 1+ J/ 2 ] = y ^^£ u (a,,p)A(z,x,k)^ drexp{ir[ W -(k-p/2)-p]} 

/dpdto f°° 
(2ir) d+1 ^ 22 ^' P) A ( Z ' x ' k ) y ^exp{ir[w- (k + p/2) • p]} 

— r 3 ^ r .Ri 1 (w,p)A(z,x,k) y^ dr exp{zr[cu - (k - p/2) • p]} 

^22 (ttf, p)A(z, x, k) y drexp{ir[u;- (k-p/2) -p]} 

/• fin((k - p/2) • p, p) + fi 22 ((k - p/2) • p, p) x dp 
7 2 lZ ' X ' J (27r) d 

+ | ( ^/n^P)-^(^P) A (, ;X;k) |^ r expMa,-(k- P /2).p]}sgn(r) = A + 5 

with 

y jWgg^k - P ) + iM^, k - P ) dp 

A -y 2 A(z ' x ' kj (2^p 

IS 



and 



B 



dpduj Rn(u, p) - R22 (w, p) 
(27r) d + 1 2 

dp -Rn(r,p) - -R 2 2(^,P 



drexp{ir[w — (k — p/2) • p]}sgn(r)A(,z, x, k) 
exp{— ir(k — p/2) • p}sgn(r)A(2i, x, k). 



Hence (j89|) indeed holds and the proof of Lemma lM. 21 is complete. 

B The proof of Lemma 13.31 

The proof is very similar to what is presented in [§] . We highlight the main differences here and refer 



the reader to that work for additional details. We let fi(z, X, K) 6 <S(] 



p2d TU>2d\ 



be a test function 



independent of V\ t 2, where X = (xi,X2), and K = (ki,k2). We define an approximation 

H £ (z, X, K) = fi(z, X, K) + X, X/e, K) + efi 2 (z, X, X/e, K). 

We will use the notation y\{z, X, K) = fj,x(z, X, X/e, K) and X, K) = [12(2, X, X/e, K). The 
functions [i\ and /i2 are to be determined. We now use (|4U|> to get 



h=0 



^ W ® W, + kJ • V £> j ^ ( 91 ) 



/' 



3=1 

\ j=i 



3=1 



s 



3=1 



^ 2 + ( ^ + ^ ^ ' VxJ I + 



3=1 



The above expression is evaluated at £ ■ = Xj/e. The term of order e 1 in L> £ vanishes since [i is 
independent of V and the fast variable £. We cancel the term of order the same way as in 

the proof of Lemma lH. 21 bv defining fMi as the unique mean-zero (in the variables V and £ = (£i, £2)) 
solution of 



/i = 0. 



(92) 



3=1 



3=1 
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It is given explicitly by 



A*(X,€,K,V) = - Hdre^ { [ §^e^^^ - f ,k 2 ) 
1 Jo JR d \ 2 ^) 2 



^) e -( kl -p) + ^ 1 .p) /i(kl + P ka) 



+ 



^l(P) e ir(k 2 .p) + ^ 2 .p) (k k 



dre r< 

^ e -( k2 -P) + ^-P) Mkl)k2 + |) 
Let us also define /i 2 as the mean zero with respect to 7iy solution of 

2 



/'i 



Mi, 



(93) 



where f = j diryf- 

In order to finish the proof of Lemma 13.31 we have to compute 



= -E 
and verify that 



/'i 



} = h + I 2 (94) 

||£* £ -£*0£*|| i2 ^- 2 ^O (95) 

as e — > 0. This is done by a straightforward but tedious calculation. We present some of the details 
for the convenience of the reader. The first term in (|94j) is 



* [Jul* ( 27 2 Jm 

Now we further split 



§^e^V l( k 1 + f,k 2 ) \=I 11+ h, 

R d (2TT) a 2 



hi = hm + hin + hni + A122 

according to the four terms in the expression for /x 2 . We compute the first and the third terms as 
they illustrate the general picture: 



I 



mi 



dV M e i(P^)i I dre rC 



i \J Rd (2ir) d 
dr [ J R 11 (^,p)e 4r ^-( kl - p / 2 > P V(ki,k 2 



^W e ir((k 1 -p/2).q)+i« 1 .q) Ai ( kl _ E _ S,^) 

(2ir) d PV 2 2 1 



(2vr) c 



(96) 



and 



'1112 



-E 



M. d (2vr) d 



e i(p-€i) / dre 



dFi(q) 



^°°dr | i? 11 ( w ,p) e l P^^)/V^-( k2 -P))/x(k 1 -|,k 2 + | 



' e ir(k 3 .q)+i« a .q) , k P k 

(27r) d PV 2 2 y 

P p, dpdu 



2 ' (2vr) 



(97) 
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The terms as in (|9*fi|l combine exactly to be equal to C* (g> C* The terms as in (|9~7|l vanish as e — > 
in the L 2 -sense - this is verified as in [B]. Notice that the a priori regularity of the Wigner measure 
in L 2 (M. 2d ) resulting from the mixture of states is crucial to obtain convergence to in (|97j) : see the 
difference between (Sj and [B] . This completes the sketch of the proof of Lemma 13.31 
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